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ABSTRACT 

We present photometric and spectroscopic data of the type II-P supernova 2003gd, 
which was discovered in M74 close to the end of its plateau phase. SN 2003gd is the first 
type II supernova to have a directly confirmed red supergiant progenitor. We compare 
SN 2003gd with SN 1999em, a similar type II-P supernova, and estimate an explosion 
date of 18th March 2003. We determine a reddening towards the supernova of E(B — 
V) = 0.14 ±0.06, using three different methods. We also calculate three new distances 
to M74 of 9.6 ± 2.8 Mpc, 7.7 ±1.7 Mpc and 9.6 ± 2.2 Mpc. The former was estimated 
using the Standardised Candle Method (SCM), for type II supernovae, and the latter 
two using the Brightest Supergiants Method (BSM). When combined with existing 
kinematic and BSM distance estimates, we derive a mean value of 9.3 ± 1.8 Mpc. 
SN 2003gd was found to have a lower tail luminosity compared to other "normal" type 
II-P SNe bringing into question the nature of this supernova. We present a discussion 
concluding that this is a "normal" type II-P supernova which is consistent with the 
observed progenitor mass of 8+2 M Q . 

Key words: stars: evolution - supernovae: general - supernovae: individual: 
SN 2003gd - galaxies: individual: M74 - galaxies: distances and redshifts 



1 INTRODUCTION 

Supernova (SN) 2003gd was discovered by Rev. Bob Evans 
on June 12.82 2003 UT (JD 2452803.32) with an apparent 
magnitude of ~13.2 mag, situated on the southern spiral 
arm of M74. The SN was soon confirmed by R. H. Mc- 
Naught on June 13.84 UT who gave a precise position of 
R. A. = l' i 36 m 42 s .65, Dec = +15°44'20".9 (Evans & Mc- 
Naught 2003). Garnavich & Bass (2003) obtained a near- 
infrared spectrum of the SN on June 13.46 UT showing 
broad Paa and Pa7 which are indicative of a supernova of 
type II. Kotak et al. (2003) confirmed 2003gd as a type II 



with an optical spectrum and estimated that it was roughly 
2 months post-explosion. Owing to the relative proxim- 
ity of M74 (9.3 Mpc; see §4.3) and its favourable inclina- 
tion, a monitoring programme was initiated. M74 was also 
the host galaxy of the recent type Ic SN 2002ap (Mazzali 
et al. 2002; Smartt et al. 2002b; Gal- Yam et al. 2002). Van 
Dyk et al. (2003) presented early BVRI light curves and 
colour curves of SN 2003gd and compared them with those 
of SN 1999em. They photometrically classified 2003gd as 
a type II-P SN and also gave estimates of the explosion 
date, 17th March 2003, and of the line of sight reddening, 
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E(B — V) — 0.13 ± 0.03. Using ground based astrometric 
calibrations Smartt et al. (2003b) and Van Dyk et al. (2003) 
estimated the position of SN 2003gd on pre-discovery Hub- 
ble Space Telescope (HST) images of M74, suggesting possi- 
ble candidates for the progenitor but both studies required 
follow-up HST images to determine the SN position with 
sufficient precision to unambiguously identify any precursor 
star. 

After observing SN 2003gd with the Advanced Camera 
for Surveys (ACS) as part of our ongoing HST programme 
on supernova progenitors, we presented the discovery of a 
star which exploded as SN 2003gd (Smartt et al. 2004). The 
progenitor star was also identified on pre-explosion Gemini 
images and the combined VRI magnitudes were used to es- 
timate its intrinsic colour and luminosity. These were found 
to be consistent with those of a red supergiant, and compari- 
son with stellar evolutionary tracks suggested that it had an 
initial main sequence mass of 8I2 Mq . This is only the third 
time that a progenitor of a genuinely confirmed SN has been 
identified and is the first apparently normal red supergiant 
to be associated with a SN of type II-P. The determination 
of intrinsic colour and luminosity require the reddening and 
distance to be known. We summarised the results of our 
measurements for both of these quantities in Smartt et al. 
(2004) and in this paper we present the full photometric and 
spectroscopic data set. Given the discovery of a red super- 
giant progenitor star, which have long been assumed to be 
the progenitors of SNe II-P, it is essential that the SN itself 
is studied and quantified. 

In this paper we present extensive photometric and 
spectroscopic data of SN 2003gd in §2 followed by an anal- 
ysis of the photometry in §3.1, where an explosion date is 
estimated. We estimate the reddening towards the super- 
nova, using three different methods, in §3.2 and analyse the 
velocity evolution in §3.3. The distance to M74 is not well 
known and as yet no Cepheid distance exists. We attempt to 
improve the situation in §4 by estimating the distance using 
two different methods and compiling all distances within 
the literature. Using the distance found we then calculate 
the amount of 56 Ni synthesised in the explosion in §5. In 
§6 the observed properties of SN 2003gd are compared to 
the semi-analytical model of Zampieri et al. (2003) and in 
§7 we discuss the nature of the supernova and the implica- 
tion for the progenitor. We then conclude in §8. Throughout 
this work we have assumed the galactic reddening laws of 
Cardelli, Clayton & Mathis (1989) with R v = 3.1. 



2 OBSERVATIONS 

2.1 Photometry of SN 2003gd 

BVRI photometry was obtained shortly after discovery and 
covers a range of ~92-169 days after explosion with a few 
late-time epochs around 490 days. These data were collected 
from several telescopes: the 1.0 m Jacobus Kapteyn Tele- 
scope (JKT), the 2.5 m Isaac Newton Telescope (INT) and 
the 3.5 m Telescopio Nazionale Galileo (TNG), all on La 
Palma, and the 1.82 m Copernico Telescope (CTA) at Asi- 
ago, Italy. The latest data were taken with Unit Telescope 
1 (UT1), Antu, of the ESO Very Large Telescope (VLT), 
Paranal, Chile. A summary of these observations can be 



found in Table 1 as well as the results from the optical pho- 
tometry. 

The JKT observations where taken with its 2048 x 2048 
SITel JAG-CCD camera and were reduced using standard 
techniques within iraf. The frames were debiased and flat- 
fielded using twilight sky flats generally taken on the night 
of the observations. In some cases flats for our required fil- 
ters were not available from that night and we employed the 
best set of flats closest in time to the supernova observations. 
An extensive photometric sequence of stars around M74 has 
been calibrated by Henden (2002) originally intended for 
monitoring of SN 2002ap. On the JKT images there were 
normally nine stars suitable for calibration of the supernova 
magnitudes (Henden numbers 100, 109, 149, 153, 156, 158, 
160, 161) and these were used to determine the zero-points 
and transformation coefficients on each night. In some cases 
the telescope pointing was slightly different or the brightest 
stars in the sequence were saturated due to increased ex- 
posure time for the fading supernova and hence not all the 
standards could be used. A minimum of three standard stars 
were always available on each night. The BVRI magnitudes 
were obtained using aperture photometry and the errors are 
the standard deviations of the individual values of the su- 
pernova magnitudes determined from each calibration star. 

The INT observations were taken with the Wide Field 
Camera (WFC) at the prime focus. The WFC consists of 
4 thinned EEV 2kx4k CCDs and SN 2003gd and its cali- 
bration stars were placed on CCD4 of the mosaic. The data 
were reduced using the standard pipeline processing applied 
to all WFC images in the Cambridge archive (Irwin & Lewis 
2001). Images were taken with the Harris V filter and the 
SDSS-like filters r' and i' . We attempted to transform the 
r' and i' aperture magnitudes to the RI system of Henden 
using 15 stars from the Henden catalogue. We derived good 
fits between the instrumental magnitudes and the intrinsic 
(R — I) and (V — R) colours for the stars and tested these 
on four other stars in common to the Henden and WFC 
catalogues. The transformations produced VRI magnitudes 
which were accurate to a mean of 0.03. However when we 
transformed the supernova CCD r'i' magnitudes to the RI 
system we got very unsatisfactory results. On 2003 August 
25, we derived R = 16.20 and I = 16.34, which are very 
discrepant from our AFOSC magnitudes taken the previ- 
ous night. As the transformations for field stars proved to 
be reliable the only explanation for this serious discrepancy 
is the strong emission line SN spectrum producing system- 
atic differences in the colour terms compared to the stellar 
SEDs. Hence we quote the V band magnitude only for these 
two nights as the r'i' transformations cannot be relied upon 
to produce magnitudes consistent with our previous calibra- 
tions. The errors in the magnitudes were calculated as with 
the JKT data. 

The CTA observations were taken with the Asiago Faint 
Object Spectrograph and Camera (AFOSC), a focal reducer 
instrument, which is equipped with a Tektronix 1024x1024 
CCD with a pixel sampling of 0.47 arcsec pixel -1 . The TNG 
observations were instead secured with the Device Opti- 
mized for LOw RESolution (DOLORES), which is equipped 
with a Loral 2kx2k CCD with a pixel scale of 0.275 arcsec 
pixel -1 . The data were reduced using standard techniques 
within iraf. The SN magnitudes were measured using the 
iraf point-spread-function (PSF) fitting task daophot. 
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Figure 1. BVRI light curves of SN 2003gd which have been 
arbitrarily shifted in magnitude for clarity. 



This procedure allows the simultaneous fitting and subtrac- 
tion of the galaxy background. For cases in which the seeing 
is fair, where the SN is relatively bright and its PSF well- 
sampled, it has been found that this method produces re- 
sults in excellent agreement with the template subtraction 
method (see Rigon et al. 2003). The errors in the derived 
magnitudes were calculated by combining in quadrature the 
standard deviation of the magnitudes found from each of 
the calibration stars, the PSF fitting error and the standard 
deviation of the magnitudes of artificial stars placed around 
the SN position having the SN magnitude. The photometry 
results, from all telescopes, are tabulated in Table 1. The 
earlier data are plotted in the BVRI light curves in Fig- 
ure 1 whereas the late-time data are shown in Figure 11. 



2.2 Spectroscopy of SN 2003gd 

Optical spectroscopic observations were taken at five tele- 
scopes and cover eight epochs, which we estimate are be- 
tween 87 and 157 days after explosion, and one late-time 
epoch at 493 days. The log of the observations is shown 
in Table 2. The data were reduced using standard routines 
within IRAF. The frames were debiased, flat-fielded and ex- 
tracted and then wavelength calibrated using Cu-Ar and Cu- 
Ne lamp spectra. The wavelength calibration was checked by 
determining the positions of the night sky lines and small 
adjustments were made. The spectra were then flux cal- 
ibrated using spectrophotometric flux standards observed 
with the same instrumental setup. The slit-widths employed 
were generally between 1 — 2.1 arcsec, hence the flux cali- 
bration is unlikely to give an accurate absolute scale. We 
used the BVRI photometry to adjust the absolute scale of 
the spectra at each epoch. Spectral BVRI magnitudes were 
calculated from the spectra using the iraf package SYN- 
phot within STSDAS. These magnitudes were compared to 
the observed photometric magnitudes and a scaling factor 
was calculated. The spectra were then read into the spectral 
analysis program DIPSO (Howarth, Murray, Mills & Berry 
2003) for further analysis. The earlier spectral evolution of 
SN 2003gd is shown in Figure 2 and is compared to that 
of the well observed SN 1999em in Figure 3. The late-time 



spectrum is shown in Figure 13 of §3.4. These spectra will 
be available through the SUSPECT 1 website. 

2.3 HST ACS and WFPC2 Observations of M74 

2.3.1 A CS HRC Imaging of SN 2003gd 

BV and I imaging of SN 2003gd was acquired with the High 
Resolution Camera (HRC) of the Advanced Camera for Sur- 
veys (ACS) on HST. SN 2003gd was observed at an epoch 
~137d post explosion as part of program G09733. These 
observations were acquired in order to image the surround- 
ings of the supernova to determine a value of the reddening 
towards stars near the supernova and to determine the ex- 
act location of the supernova with respect to nearby stars. 
These observations allowed the location of the supernova to 
be determined on pre-explosion observations and a direct 
identification of the progenitor star (Smartt et al. 2004). 
Details of the HST ACS observations are given in Table 3. 

The on-the-fiy re-calibrated (OTFR) ACS images were 
obtained from the Space Telescope European Coordinating 
Facility archive. Photometry was conducted on these frames 
using the iraf package daophot (Stetson 1987) and its in- 
corporated PSF fitting algorithm ALLSTAR (Stetson & Harris 
1988). HST photometric system magnitudes, of stars and the 
supernova in these images, were converted to the Johnson- 
Cousins' system using the transformations of Holtzman et al. 
(1995). 

2.3.2 WFPC2 Observations of M74 

M74 was imaged, in B and J, with the Wide Field Plan- 
etary Camera 2 (WFPC2) instrument on HST as part of 
program GO9042 although these observations did not in- 
clude the location of SN 2003gd. Details of these obser- 
vations are given in Table 3. The OTFR images were re- 
trieved from the Space Telescope European Coordinating 
Facility archive, with cosmic ray split frames pre-combined. 
Photometry was conducted on these frames using the HST- 
phot version 1.15b package (Dolphin 2000a, b). HSTphot in- 
cludes corrections for charge transfer efficiency, PSF vari- 
ation and aperture size, and converts instrumental magni- 
tudes to standard UBVRI magnitudes automatically. Stars, 
suitable for the brightest supergiants distance determina- 
tion technique (see §4.2), were selected from the photom- 
etry output using the HSTphot object type classification 
scheme. Objects which were classified by HSTphot as ex- 
tended, blended or containing bad pixels were discarded, as 
were probably blended stars or stars in crowded areas with 
PSF fit x 2 > 2.5. 



3 ANALYSIS 

3.1 Explosion epoch 

The first photometric data were obtained ~6 days after dis- 
covery. As can be seen from the BVRI light curves (Figure 
1) SN 2003gd was discovered close to the end of its plateau 

1 http://bruford.nhn.ou.edu/~suspect/ 
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Table 1. Journal and results of optical photometry of SN 2003gd 



Date 




JD 

(2450000+) 


Phase 
(days) 


B 






V 






R 




Telescope + 
Instrument 


Observer 


2003 


Jun 19 


2809 


20 


92 


15.47(0.02) 


14. 


.21(0. 


.02) 


13 


.68(0 


,02) 


13.34(0.03) 


JKT+JAG 


Jonker 


2003 


Jun 20 


2810 


21 


93 


15.50(0.02) 


14, 


.18(0. 


.04) 


13 


,63(0 


,02) 


13.30(0.04) 


JKT+JAG 


Jonker 


2003 


Jun 21 


2811 


21 


94 


15.52(0.02) 


14 


.24(0. 


.04) 


13 


,73(0 


,05) 


13.33(0.03) 


JKT+JAG 


Jonker 


2003 


Jun 22 


2812 


.21 


95 


15.53(0.01) 


14 


.23(0. 


.01) 


13 


,69(0 


,02) 


13.38(0.05) 


JKT+JAG 


Jonker 


2003 


Jun 23 


2813 


20 


96 


15.57(0.02) 


14, 


.24(0. 


.01) 


13 


.71(0 


,01) 


13.35(0.02) 


JKT+JAG 


Jonker 


2003 


Jun 27° 


2817. 


.66 


100 


15.58(0.03) 


14 


.29(0. 


.03) 


13 


,75(0. 


,02) 


13.44(0.04) 


TNG+DOLORES 


Zampieri 


2003 


Jun 28 




ly 


101 


15.62(0.01) 


14 


.37(0. 


.04) 


13 


,76(0 


,02) 


13.35(0.05) 


JKT+JAG 


Vermaas 


2003 


Jul 02 


2822 


20 


105 


15.75(0.02) 


14, 


.42(0. 


.07) 


13 


.82(0 


,04) 


- 


JKT+JAG 


Vermaas 


2003 


Jul 03 


ZoZo 


i o 

iy 


106 


15.78(0.01) 


14, 


.45(0. 


.06) 


13 


.85(0 


.04) 


- 


JKT+JAG 


Vermaas 


2003 


Jul 08 


2828. 


.18 


111 


15.98(0.06) 


14, 


.59(0. 


,03) 


13 


,99(0 


.08) 


13.63(0.03) 


JKT+JAG 


Exter 


2003 


Jul 11 


2831 


21 


114 


16.15(0.02) 


14, 


.78(0. 


,03) 


14 


,14(0 


.04) 


13.76(0.05) 


JKT+JAG 


Extcr 


2003 


Jul 12 


ZooZ. 


zz 


115 


16.24(0.02) 


14 


.84(0 


,03) 


14 


.21(0 


,04) 


13.82(0.04) 


JKT+JAG 


Exter 


2003 


Jul 14 


zoo4 


.zz 


117 


16.46(0.02) 


14, 


,98(0, 


,03) 


14 


,35(0 


.04) 


13.93(0.06) 


JKT+JAG 


Extcr 


2003 


Jul 15 


Zooo 


lo 


118 


16.58(0.04) 


15, 


,04(0, 


,05) 






- 


- 


JKT+JAG 


Exter 


2003 


Jul 16 


ZOOO 




119 


16.66(0.05) 


15, 


,17(0 


.06) 


14 


,56(0 


,06) 


14.11(0.09) 


JKT+JAG 


Extcr 


2003 


Jul 19 


2839 


23 


122 


17.35(0.10) 


15, 


,54(0, 


.06) 






- 


- 


JKT+JAG 


Exter 


2003 


Jul 21 


OQ/l 1 

zo41. 


OD 


124 


17.72(0.04) 


16 


,02(0, 


,04) 


15 


,23(0 


,03) 


14.59(0.05) 


CTA+AFOSC 


Pastorello 


2003 


Jul 23 


2843. 


18 


126 


- 


16. 


,75(0, 


,04) 


15 


,77(0 


,04) 


- 


JKT+JAG 


Verley/Leon 


2003 


Jul 24 


2844. 


19 


127 


- 


16. 


,89(0. 


,03) 


16 


,03(0 


,03) 


- 


JKT+JAG 


Verley/Leon 


2003 


Jul 26 


2847 


19 


128 


- 


17. 


.24(0. 


02) 


16 


,30(0 


,02) 


- 


JKT+JAG 


Verley/Leon 


2003 


Jul 27 


2848 


18 


129 


- 


17. 


.26(0. 


01) 


16 


,34(0 


,02) 


- 


JKT+JAG 


Verley/Leon 


2003 


Jul 28 


2849. 


.18 


130 


- 


17. 


.27(0. 


.01) 


16 


,34(0 


,03) 


- 


JKT+JAG 


Verley/Leon 


2003 


Jul 29 


2850. 


18 


131 


- 


17. 


.36(0. 


.02) 


16 


,38(0 


,02) 


- 


JKT+JAG 


Verley/Leon 


2003 


Jul 30 


2851. 


.18 


132 


- 


17. 


.37(0. 


.04) 


16 


,39(0 


,02) 


- 


JKT+JAG 


Verley/Leon 


ZUUo 


Tnl "il 
j ui or 


2852 


18 


loo 




1 7 


17ffi 




ID. 


iQ(n 
oy \yj 


UO ) 




TT\ T_l_ TAP 


Verley/Leon 


2003 


Aug 07 


2858 


73 


141 




17. 


44(0. 


03) 








15.77(0.04) 


TNG+DOLORES 


Zampieri 


2003 


Aug 08 b 


2859. 


72 


142 


19.24(0.05) 








16 


,36(0 


,03) 




TNG+DOLORES 


Zampieri 


2003 


Aug 24 


2875 


51 


158 


19.34(0.10) 


17. 


,62(0, 


.05) 


16 


,75(0 


,05) 


16.16(0.05) 


CTA+AFOSC 


Bcnctti 


2003 


Aug 25 


2876 


.71 


159 




17. 


,56(0. 


05) 










INT+WFC 


Irwin 


2003 


Sep 04 


2886. 


.61 


169 




17. 


.64(0. 


,05) 










INT+WFC 


Irwin 


2004 


Jul 15 


3201 


69 


485 










20 


.55(0 


,09) 




TNG+DOLORES 


Zampieri 


2004 


Jul 23 


3209 


65 


493 




21. 


.33(0. 


.26) 








20.30(0.24) 


TNG+DOLORES 


Zampieri 


2004 


Jul 23 


3209 


.91 


493 


21.76(0.06) 


21 


.32(0. 


.07) 


20 


,59(0 


,09) 


20.31(0.04) 


VLT-UT1+FORS2 


Patat/Pignata 



NOTE:- Figures in brackets give the errors associated with the magnitudes. A full account of these errors can be found in §2.1. 
a U = 17.19, b U = 21.39 

JKT — 1.0 m Jacobus Kaptcyn Telescope, La Palma 
TNG — 3.5 m Tclcscopio Nazionalc Galileo, La Palma 
CTA — 1.82 m Copcrnico Telescope, Asiago, Italy 
INT — 2.5 m Isaac Newton Telescope, La Palma 

VLT-UT1 = 8.2 m ESO Very Large Telescope, Antu (UT1J, Paranal. Chile 

Table 2. Journal of spectroscopic observations of SN 2003gd 



Date 


JD 


Phase 


Range 


Resolution 


Instrument 


Telescope 


Observer 




(2450000+) 


(days) 


(A) 


(A) 








2003 Jun 14 


2804.72 


87 


3953-7499 


4 


ISIS 


WHT 


Benn 


2003 Jun 25 


2815.68 


98 


3402-9998 


8 


IDS 


INT 


Lcnnon 


2003 Jun 27 


2817.67 


100 


3500-10430 


~14 


DOLORES 


TNG 


Zampieri 


2003 Jul 09 


2829.68 


112 


3200-11000 


8 


IDS 


INT 


Prada 


2003 Jul 21 


2841.60 


124 


3494-7809 


~24 


AFOSC 


CTA 


Pastorello 


2003 Aug 07 


2858.71 


141 


3156-8066 


~14 


DOLORES 


TNG 


Zampieri 


2003 Aug 22 


2873.72 


155 


3147-10440 


~14 


DOLORES 


TNG 


Zampieri 


2003 Aug 24 


2875.55 


157 


3486-7475 


~24 


AFOSC 


CTA 


Benetti 


2004 Jul 23 


3209.90 


493 


4233-9632 


~13 


FORS2 


VLT-UT1 


Patat/Pignata 



WHT = 4.2 m William Hcrschcl Telescope, La Palma 
TNG — 3.5 m Tclcscopio Nazionale Galileo, La Palma 
CTA — 1.82 m Copcrnico Telescope, Asiago, Italy 
INT — 2.5 m Isaac Newton Telescope, La Palma 

VLT-UT1 = 8.2 m ESO Very Large Telescope, Antu (UT1J, Paranal, Chile 
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Figure 2. Spectral evolution of SN 2003gd from 87-157 days after explosion, where the spectra have been arbitrarily shifted in flux for 
clarity. 



Table 3. Summary of HST observations of SN 2003gd and M74. 



Date 


Filter 


Dataset 


Exposure (s) 


Instrument 


2003 Aug 1 


F439W 


J8NV01011/21 


2500 


ACS/HRC 


2003 Aug 2 


F555W 


J8NV01031/41 


1100 


ACS/HRC 


2003 Aug 2 


F814W 


J8NV01051 


1350 


ACS/HRC 


2001 Oct 3 


F450W 


U6EA0101/02 


460 


WFPC2 


2001 Oct 3 


F814W 


U6EA0103/04 


460 


WFPC2 



phase when the host galaxy became visible after its con- 
junction with the sun. The plateau phase ends ~30 days 
later when the nebular tail phase begins. The light curve of 
SN 2003gd was compared to that of SN 1999em, a similar 
type II-P supernova, using a x 2 -htting algorithm to adjust 
the time and apparent magnitude to find the best fit. The 
SN 2003gd data from after JD 2452845 were not included in 
the fit due to the difference in tail luminosities. The BVRI 
best fits are shown in Figure 4 with the shift in time and 
apparent magnitude inset in each figure. The results from 
the x 2 -fitting algorithm are shown in Table 4 as well as the 
weighted average of the shift in days for the VRI bands. Al- 
though B has a large reduced-^ 2 , the fit to the "knee" in the 
light curve looks as good as the others. The poor fit could be 
due to one of two reasons: i) under estimation of the errors 
in the photometry or ii) the errors being non-Gaussian due 
to the lower magnitude in the B-band. The B-band data 



was subsequently removed from the explosion date calcula- 
tion, but had it been included there would have been little 
difference in the result. The errors in At, the shift in days, 
and Am, the shift in magnitude, are the errors determined 
from the x 2_n t at the la level and are shown in brackets 
in Table 4. The explosion date was determined using the 
weighted average of At, in VRI, and the explosion epoch of 
SN 1999em (Hamuy et al. 2001), which was estimated as JD 
2451478. 8±0.5, within two days of Elmhamdi et al. (2003b). 
We find the explosion date of SN 2003gd therefore to be JD 
2452717±21, which corresponds to 18th March 2003. The 
error in the explosion epoch was estimated from the error 
in the weighted average of At, the error in the SN 1999em 
explosion date and the uncertainty in the duration of the 
plateau of SN 2003gd. This latter uncertainty was estimated 
from the observed parameters of a sample of 13 SNe II-P 
from Hamuy (2003, Table 3) where the mean plateau du- 
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99em + 79d_ 




03gd +98d 




Figure 3. Comparison of the spectra of SN 2003gd and SN 1999em at similar phases, where the spectra have been arbitrarily shifted in 
flux for clarity. The spectra of both supcrnovae show remarkable similarities at all epochs. 



Table 4. Results from the x 2 -fitting algorithm which adjusts the time, At, and apparent magnitude, Am, of the "model" light curve 
(SN 1999em) to find the best fit to the data points of SN 2003gd, where v = number of data points - number of degrees of freedom. 



Filter 


Rcduced-x 2 


V 


At (days) 


Am 


B 


6.92 


14 


1236.20 


-0.39 


V 


1.37 


12 


1239.65(0.53) 


-0.01(0.02) 


R 


1.02 


12 


1237.39(0.65) 


-0.05(0.02) 


I 


0.69 


10 


1238.11(1.14) 


-0.08(0.03) 


VRI 






1238.67(0.39) 





ration was found to be 131±21 days. The most significant 
error is introduced by the assumption that the plateau of 
SN 2003gd was of a similar duration to SN 1999em. We have 
based this assumption on the spectral similarity of the two 
events (see Figure 3 and §4.1) and their comparable masses 
(Smartt et al. 2002a; Elmhamdi et al. 2003b; Smartt et al. 
2004). The error of ±21 days amply accounts for any reason- 
able intrinsic difference that may exist. Given this explosion 
epoch we estimate that the supernova was discovered ~86 
days after explosion. The light curves of both supernovae are 
comparable during the plateau phase, but SN 2003gd is sub- 
stantially fainter than SN 1999em in the late-time nebular 
phase. This is most likely due to a lower 56 Ni mass synthe- 
sised in the explosion, which was also suggested by Van Dyk 
et al. (2003). 



3.2 Reddening Estimation towards SN 2003gd 

The reddening was estimated towards SN 2003gd using three 
different methods. It was first of all estimated using the 
colours of the supernova compared to SN 1999em; secondly 
using HST photometry of the surrounding stars and finally 
from the reddening towards nearby Hn regions. 



3.2.1 Comparison of the colour evolution of SN 2003gd 

This method assumes that SNe II-P all reach the same in- 
trinsic colour towards the end of the plateau phase. This 
is, in turn, based on the assumption that the opacity of 
SNe II-P is dominated by electron scattering and, there- 
fore, the supernovae should reach the temperature of hy- 
drogen recombination at this point (Eastman et al. 1996; 
Hamuy 2004a). As shown in §2.2 the two supernovae are 
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Figure 4. BVRI light curves of SN 2003gd over-plotted with the light curve for SN 1999cm from Hamuy ct al. (2001). Filled circles are 
from data in Table 1 and open circles are from Van Dyk et al. (2003) for comparison. The reduced-x 2 of the fit for BVRI are as follows: 
6.92, 1.37, 1.02 and 0.69. 



strikingly similar at all epochs of our sample. It is therefore 
reasonable to assume that SN 2003gd would indeed reach the 
same colour as SN 1999em. Hamuy (2004a) investigated this 
method using a sample of 24 SNe II-P with BVI photome- 
try using SN 1999em as the comparison supernova. Hamuy 
(2004a) found serious discrepancies between the results ob- 
tained from (B — V) and (V — I) colours and notes that this 
is an unsatisfactory technique and other methods should be 
explored in future. 

Bearing this in mind we carried out the analysis us- 
ing (B — V), (V — R) and (V — I). These colour curves for 
SN 2003gd can be seen in Figure 5. SN 1999em was exten- 
sively modelled by Baron et al. (2000) who determined a 
hard upper limit of E(B — V) < 0.15 and estimated a likely 
reddening of E(B—V) ~ 0.05—0.10. In this paper we adopt a 
reddening of E(B - V) = 0.075 ± 0.025 for SN 1999em. The 
colour curves of SN 1999em were firstly dereddened using 
this value and were then shifted in time using the weighted 
average of At discussed in §3.1. A x 2 -fitting algorithm was 
then used to compare the colours of SN 2003gd with those 
of SN 1999em for JD < 2452830, the end of the plateau 
phase. The results of the x 2 ~fit are listed in Table 5 and 
shown in Figure 6. The errors in the colour excess, shown 
in brackets, were estimated from the x 2_n t at the la level. 
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Figure 5. Observed (B — V), (V-R), (V - 1) and (fl-J) colour 
curves of SN 2003gd plotted against the time since explosion. 
Open circles show direct observed colours whereas closed symbols 
are colours determined from the interpolated light curve. 



The systematic error, that the uncertainty in the reddening 
of SN 1999em introduces, was excluded from the fit to be 
added later. 
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Table 5. Results of x 2 -fitting algorithm to determine the reddening of SN 2003gd, where v = number of data points - number of degrees 
of freedom. 



Colour 


Rcduced-x 2 


V 


Colour Excess 


A V 


E(B - V) 


B-V 


2.20 


9 


-0.16 


-0.49 


-0.16 


V-R 


0.27 


9 


0.01 


0.06 


0.02 


V -I 


1.01 


7 


0.13(0.02) 


0.31(0.05) 


0.10(0.02) 



LO 
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: A v = -0.49 






- A v = 0.06 


i - 






1 1 1 1 1 1 1 

- A y = 0.31 


H 1 1 1 1 

I 



2750 2800 2850 



JD(2450000+) 

Figure 6. (B - V), (V - R) and (V - I) colour evolutions 
of SN 2003gd compared with the intrinsic colour evolutions of 
SN 1999em (Hamuy ct al. 2001). The visual extinction required 
for each fit are given in the panels. The reduced-x 2 for (B — V), 
(V - R) and (V - I) arc 2.20, 0.27 and 1.01, respectively. Open 
circles show direct observed colours whereas closed symbols are 
colours determined from the interpolated light curve. 



Hamuy (2004a) finds the (V — I) reddening estimate to 
be better behaved than the (B — V) estimate, and in the case 
of SN 2003gd we find this also to be true. The (B - V) fit 
yields a negative reddening which is unphysical and could be 
due to line blanketing as suggested by Hamuy. The (V — R) 
estimate also gives an unrealistically small reddening value, 
which is less than the Galactic reddening of E(B — V) = 0.07 
(Schlegel et al. 1998). In any case the reduced-x 2 values for 
the (B — V) and ( V — R) are indicative of poor fits and were 
excluded from the reddening estimate. Using the (V — 7) 
colour fit we estimate the reddening to be E(B — V) = 
0.10 ±0.03, which is comparable to that determined by Van 



-.5 



-.4 - 




'-.4 -.3 -.2 -.1 .1 .2 .3 

B-V 

Figure 7. A two colour diagram showing 25 stars, of M74, within 
6 arcsec of the location of SN 2003gd. The stars have been dered- 
dened by E(B — V) = 0.13. The intrinsic colour sequence of su- 
pcrgiants (Drilling & Landolt 2000) is shown as the solid line. 
The reddening vector is indicated by the arrow. 

Dyk et al. (2003). The error was determined by combining 
in quadrature the error from the x 2-n t and the systematic 
error in the reddening of SN 1999em. The error that the 
uncertainty in the explosion date introduces was found to 
be negligible. 

3.2.2 Reddening towards the neighbouring stars 

Three colour ACS photometry (see §2.3.1) was used to es- 
timate the reddening towards SN 2003gd. (B — V) and 
(V — I) colours of 25 stars within 6 arcsec of SN 2003gd were 
compared with the intrinsic supergiant colour sequence of 
Drilling & Landolt (2000). The reddening was calculated us- 
ing a x 2 -minimisation of the displacement of the stars from 
the intrinsic supergiant colour sequence, for a range of values 
of E(B - V). The reddening vector, in the (B - V)/(V - I) 
colour plane, assumed the reddening laws of Cardelli et al. 
(1989) with Rv = 3.1. Using this method the reddening was 
estimated as E(B - V) = 0.13 ± 0.07. The positions of the 
25 dereddened stars and the intrinsic supergiant colour se- 
quence, on the (B — V)/(V — I) colour plane, are shown in 
Figure 7. 

3.2.3 Reddening towards nearby Hll regions 

The abundance gradient in M74 has been studied by Bel- 
ley & Roy (1992), who observed the nebular lines of 132 Hll 
regions. For each Hll region the authors calculated an extinc- 
tion at H/3, which gives the total extinction i. e. a combined 
Galactic plus extragalactic extinction. The closest Hll re- 
gions to SN 2003gd, which were analysed, are Belley & Roy 
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Table 6. Nearby Hn regions to SN 2003gd in M74 analysed by 
Belley & Roy (1992). 



Hn region 


c(H/3) 


A v 


E(B - V) 


72 


0.47 


1.01 


0.33 


73 


0.06 


0.13 


0.04 


Mean 






0.19(0.15) 



Table 7. Summary of the results from the three methods used in 
determining the reddening of SN 2003gd. 



Method 


E(B-V) 


Comparison with SN 1999em 


0.10(0.03) 


HST photometry 


0.13(0.07) 


Nearby Hn regions 


0.19(0.15) 


Mean 


0.14(0.06) 



numbers 72 and 73. The logarithmic extinction c(H/3), vi- 
sual extinction Av and the colour excess E(B — V) for these 
Hn regions are given in Table 6 along with their mean. The 
visual extinction was calculated from the logarithmic extinc- 
tion and an extinction relationship of ^4861 /Av = 1.16 with 
Rv = 3.1 (Cardelli et al. 1989). 

3.2.4 Visual extinction of SN 2003gd 

A summary of the results of the methods discussed is given 
in Table 7 along with the mean of the methods, E(B — V) = 
0.14±0.06. The error on the mean was found to be unrealis- 
tically small at 0.03, therefore we have quoted the combined 
errors of the three methods. The Galactic extinction law 
with Rv = 3.1 was used to calculate the visual extinction, 
which was found to be Av = 0.43 ± 0.19. 

3.3 Expansion velocity of SN 2003gd 

We measured the expansion velocity of the ejecta of 
SN 2003gd from the minimum of the blue-shifted absorption 
trough of the Fe 11 A5169 line, as in Hamuy (2004b). This 
absorption line is asymmetric and shows signs of blending, 
possibly with Ti 11 and Mg I A5167 (Pastorello 2003). The 
absorption trough was fitted by three Gaussians using the 
Emission Line Fitting package (elf) within the STARLINK 
spectral analysis package DIPSO. The minimum was found 
from the convolution of the three Gaussians and the error 
was estimated from the difference between this and a sin- 
gle Gaussian fit. The velocities measured in this way for 
SN 2003gd are listed in Table 8. As a consistency check 
several epochs of SN 1999em, from Hamuy et al. (2001), 
Leonard et al. (2002a) and Elmhamdi et al. (2003b), were 
also measured using the same technique. The velocities mea- 
sured here are listed in Table 9. The velocity evolutions of 
SNe 2003gd and 1999em are plotted in Figure 8. 

The velocity evolution of SN 2003gd closely follows that 
of SN 1999em at all epochs. As the supernovae age it be- 
comes increasingly difficult to measure the minimum of the 
Fe 11 A5169 line as a blended trough begins appearing as 
shown in Figure 9. The velocity measurements carried out 
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Table 8. Velocities derived from the minimum of the Fe II A5169 
line for SN 2003gd. 



Date 


JD 


Phase 


VFeIIA5169 




(2450000+) 


(days) 


(km s 1 ) 


2003 Jun 14 


2804.72 


87 


2514(237) 


2003 Jun 25 


2815.68 


98 


2282(126) 


2003 Jun 27 


2817.67 


100 


2398(126) 


2003 Jul 09 


2829.68 


112 


2050(126) 


2003 Jul 21 


2841.60 


124 


1934(126) 


2003 Aug 07 


2858.71 


141 


1876(872) 


2003 Aug 22 


2873.72 


155 


1586(525) 


2003 Aug 24 


2875.55 


157 


1237(126) 



Table 9. Velocities derived from the minimum of the Fe II A5169 
line for SN 1999em using the technique described here. 



Date 


JD 


Phase 


VFcIIA5169 


Spectral 




(2450000+) 


(days) 


(km s- 1 ) 


Source 


1999 Nov 09 


1492.14 


13 


7101(271) 


1 


1999 Nov 14 


1496.67 


18 


6088(182) 


2 


1999 Nov 19 


1501.66 


23 


5268(173) 


2 


1999 Dec 15 


1527.74 


49 


3722(378) 


1 


1999 Dec 17 


1529.74 


51 


3542(294) 


1 


1999 Dec 18 


1530.70 


52 


3394(311) 


3 


1999 Dec 29 


1541.70 


63 


3250(287) 


3 


1999 Dec 31 


1543.66 


66 


3098(272) 


2 


2000 Jan 13 


1556.84 


78 


2834(189) 


1 


2000 Feb 01 


1575.74 


97 


2439(360) 


1 


2000 Feb 09 


1584.35 


106 


2173(263) 


3 


2000 Feb 16 


1590.55 


112 


2207(341) 


3 


2000 Mar 01 


1604.74 


126 


1855(569) 


1 


2000 Mar 13 


1616.51 


138 


1477(353) 


3 


2000 Mar 15 


1618.64 


140 


1306(305) 


1 



1 Leonard ct al. (2002a), 2 Hamuy et al. (2001), 3 Elmhamdi 
et al. (2003b) 
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Figure 8. Velocity evolution of SN 2003gd and SN 1999cm de- 
rived from the minimum of the blue-shifted absorption trough of 
Fe II A5169. The superscripts in the figure denote the source of the 
velocity measurements: (1) this paper, (2) Leonard et al. (2002a), 
(3) Hamuy et al. (2001). 
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Figure 9. Spectrum of SN 1999em (Leonard et al. 2002a) from 
161 days after explosion showing blend appearing in the Fe II 
A 5169 line. 
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Figure 10. Velocity evolution of SN 1999cm derived from the 
minimum of the blue-shifted absorption troughs of Fell AA4924, 
5018, 5169, taken from Leonard ct al. (2002a). 
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Figure 11. BVRI nebular and late-time photometry of 
SN 2003gd with a linear fit of the earlier nebular data, in the 
V-band, shown with a dashed line. The late-time photometry, in 
the V-band, is lower than what we would expect from the earlier 
photometry. 



100 200 300 400 500 
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Figure 12. Pseudo-bolometric (UVOIR) light curve of 
SN 2003gd with a linear fit of the earlier nebular data shown 
with a dashed line and the decay slope of 56 Co — > 56 Fe shown 
with a dotted line. 



in this work on the spectra of Hamuy et al. (2001), Leonard 
et al. (2002a) and Elmhamdi et al. (2003b) are compara- 
ble with the existing measurements. This demonstrates the 
consistency of the techniques used. Although the methods 
are consistent we propose that the Fe II A5169 line is not a 
good indication of the ejecta velocity at later epochs. Fig- 
ure 10 shows the velocities derived from the Fe II AA4924, 
5018, 5169 lines, taken from Leonard et al. (2002a), which 
are traditionally used to determine the photospheric veloc- 
ity. These authors, however, point out the consistently lower 
velocity derived from the Fe II A4924 line compared with 
those calculated from the other Fe II lines. Leonard et al. 
suggest that this is due to significant blending with Ba II 
A4934, which was also found by Pastorello (2003). Figure 
10 demonstrates this, but also shows that the Fe II A5169 
line, at around 140 days, gives significantly lower velocities. 
Although Fe II A5169 yields lower velocities at these epochs, 
for our purposes of comparison, this is of little consequence. 



3.4 Late-time photometry and spectroscopy 

Late-time photometry was obtained around 490 days after 
explosion and is plotted in Figure 11 along with the nebular 
phase photometry. The details of these observations can be 
found in Table 1. A linear fit has been applied to the early 
nebular phase data in the V-band, which is shown with a 
dashed line. The decay rate from this fit was estimated to 
be 0.98 mag 100d _1 , which corresponds to the decay rate 
of 56 Co, suggesting that little or no 7-rays escaped during 
this time. The pseudo-bolometric (UVOIR, see §7.1) light 
curve is shown in Figure 12 where the dashed line is the 
fit to the early nebular data and the dotted line is the de- 
cay slope of 56 Co — » a6 Fe. A late-time spectrum was also 
obtained, around 493 days after explosion, and is shown in 
Figure 13 alongside the penultimate spectrum of SN 2003gd 
and a spectrum of SN 1999em of a similar epoch. 

Figure 11 shows a deficit in the latest value of the V- 
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Figure 13. Late-time spectrum of SN 2003gd at 493d shown 
with its penultimate spectrum and a spectrum of SN 1999em of 
a similar epoch. 



band luminosity when compared with that which is expected 
from the linear fit. This deficit could either be due to dust 
formation or the supernova becoming more transparent to 
7-rays, allowing leakage to occur, or a combination of both. 
The fit in Figure 12, the dashed line, would suggest that the 
latest point is actually brighter than expected, but due to 
the limited data in the tail it is difficult to draw any conclu- 
sions. If we assume the 7-rays are fully thermalised, as the 
fit to the V-band data would suggest, then the luminosity 
is substantially less than what we would expect. Elmhamdi 
ct al. (2003b) noted a blue-shift of the line profiles of [OI] 
AA6300, 6364 A and Ha, and attributed this to the formation 
of dust. Figure 14 shows the profiles of Ha from three epochs; 
124d, 157d and 493d post explosion, where the vertical line 
marks the maximum of the Ha emission from the earliest 
spectrum. The peak of the Ha emission from 157d shows 
a slight red-shift, but the latest spectrum clearly shows a 
skewing towards the blue. The [OI] AA6300, 6364 A lines, 
however, are not clear enough in the earlier spectra to ver- 
ify this blue-shift. Nevertheless, the drop in luminosity at 
late-times, combined with the apparent blue-shift of the Ha 
emission, does indeed suggest dust formed between 157 and 
490 days. 




+ 157d xl.5 
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Figure 14. Ha profiles of SN 2003gd from 124d, 157d and 493d 
post explosion. The vertical line marks the maximum of the Ha 
emission from the earliest spectrum. The peak of the Ha emission 
from 157d shows a slight red-shift of the peak, but the latest 
spectrum clearly shows a shift to the blue. 



4 DISTANCE ESTIMATES TO M74 

4.1 Type II-P standard candle method distance 
estimate 

Kowal (1968) first suggested that supernovae of type la 
(SNe la) could be useful distance indicators as their peak 
magnitudes and redshifts were found to be closely corre- 
lated. From observations of high-z supernovae it is possible 
to measure the history of the expansion of the universe. In 
recent years many separate studies have been carried out 
using this SNe la method. One of the most recent studies, 
carried out by Tonry et al. (2003), confirmed the results of 
two previous studies, by Riess et al. (1999) and Perlmutter 
et al. (1999), who found that the universe is actually acceler- 
ating. There is currently great debate on whether SNe la are 
in fact standard candles and work on methods of standar- 
dising these supernovae is underway (Altavilla et al. 2004, 
and references therein). In order to verify the results from 
SNe la methods, Hamuy & Pinto (2002) investigated the 
plateau subclass of type II supernovae (SNe II-P) for use as 
standard candles. Although SNe II-P showed a wide range 
of luminosities at all epochs Hamuy & Pinto (2002) showed 
that the velocity of the ejecta of SNe II-P and their bolo- 
metric luminosities during the plateau phase are highly cor- 
related. Hamuy & Pinto (2002) formulated a Standardised 
Candle Method (SCM) for SNe II-P which can be solved for 
the Hubble Constant, provided a suitable distance calibra- 
tor is known. Hamuy in two further papers (Hamuy 2004a, b) 
firstly confirmed the SCM using a sample of 24 SNe II-P and 
then calibrated the Hubble diagram with known Cepheid 
distances to four of these supernovae. Hamuy (2004b) cal- 
culated the Hubble constant, using V and / data, to be 
H (V) = 75±7kms- 1 Mpc _1 and H Q (I) = 65 ± 12 km s" 1 
Mpc~ . In this paper we use the weighted average of these 
results, Ho = 72 ±6 km s _1 Mpc -1 , which is comparable to 
Ho = 71 ± 2 km s _1 Mpc _1 derived as part of the HST Key 
Project using SNe la (Freedman et al. 2001). This result 
demonstrates the integrity of this method. 

We use the results from Hamuy (2004b) conversely to 
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Figure 15. Velocity evolution of SN 2003gd compared with other 
similar SNe II-P (1999em, 1999gi and 1992ba) and contrasting 
SNe II-P (1999cr, 1992am and 1999br). The superscripts in the 
figure denote the source of the velocity measurements: (1) this 
paper, (2) Leonard ct al. (2002a), (3) Hamuy et al. (2001) and 
(4) Leonard et al. (2002b). 



estimate the distance to SN 2003gd, using Equations 1 and 
2: 



D(V) 



IQ^(V 50 -A V +6.249(±1. 35)Jo 9 (vbo/5000)+1. 464(±0.1B)) 



Ho 



(1) 



in^( / 50-A 7 +5.445(±0.91)(o9(, ; 5o/5000) + 1.923(±0.11)) 

D(I) = ~ 77 (2) 

Ho 

where V50, /50 and V50 are the V and I magnitudes, and the 
expansion velocity, in km s _1 , at day 50. As SN 2003gd was 
in conjunction with the sun 50 days after explosion we need 
to infer these quantities. 

The velocity evolution of SN 2003gd and SN 1999em is 
shown in Figure 8. In order to infer the velocity of SN 2003gd 
at day 50 from SN 1999em it is necessary to determine 
whether the evolution of SN 1999em is common to other 
similar type II-P supernovae or if it is in any way pecu- 
liar. Figure 15 shows the expansion velocities of SN 2003gd 
compared to other similar SNe II-P (SNe 1999em, 1999gi 
and 1992ba) and contrasting SNe II-P (1999cr, 1992am and 
1999br). In all cases, apart from SN 1999gi, the ejecta veloc- 
ities were measured from the minimum of the blue-shifted 
absorption trough of the Fe n A5169 line. In the case of 
SN 1999gi the velocities were measured from the minimum 
of the available weak, unblended line features (Leonard et al. 
2002b). SNe 1999em, 1999gi, 1992ba are all very similar 
SNe II-P with comparable velocities at day 50, V50, and to- 
wards the end of the plateau phase. The close correspon- 
dence of the velocities of these SNe with 2003gd later in the 
plateau phase (~100 days) suggests that SN 2003gd would 
also have had a similar V50. To illustrate the difference be- 
tween these and other SNe II-P, with quite contrasting ve- 
locity evolutions, we plot SNe 1999cr, 1992am and 1999br. 
Within the sample of SNe II-P studied there is no known 
example of a type II-P supernova which has a similar ve- 
locity evolution to that of SN 2003gd between 80-130 days 
which then diverges from the 1999em/1999gi/1992ba slope 
at day 50. This assumption allows us to estimate the veloc- 
ity of SN 2003gd at day 50. Due to the similarities between 
SN 2003gd and SN 1999em we take the weighted average of 



three independent measurements of W50 for SN 1999em and 
adopt this for SN 2003gd. The independent measurements 
are from Hamuy et al. (2001), Leonard et al. (2002a) and 
this paper, giving an adopted value of 1150 = 3694 ± 981 km 
s _1 for SN 2003gd, where the error is the standard deviation 
of v 50 in the sample of 22 SNe II-P from Hamuy (2004b). 
The two extreme SNe, 1992am and 1999br, were excluded in 
this standard deviation calculation due to their peculiarities. 

The apparent magnitudes for SN 2003gd at day 50 were 
estimated using Am from the x 2 -htting results (Table 4) and 
V50 and /50 magnitudes for SN 1999em (Hamuy 2004b), giv- 
ing V50 = 13.97±0.50 and J 50 = 13.27 ±0.52 for SN 2003gd. 
The errors adopted are the following errors combined in 
quadrature: 

(i) The error in the magnitude at day 50 for SN 1999em. 

(ii) The error in the X 2_n t- 

(iii) An error associated with the assumption that the 
light curves of SNe 2003gd and 1999em are intrinsically 
the same. We have estimated this from the standard de- 
viation in the absolute magnitudes of SNe 1992ba, 1999cm 
and 1999gi (Hamuy 2003), giving (M v ) = -15.92 ± 0.49, 
(Mi) = -16.35 ±0.49. 

(iv) An error associated with the plateau length, which 
was taken to be 0.09 mag in V and 0.17 mag in I. 

We use the reddening estimated in §3.2 to give Ay = 
0.43 ± 0.19 and Ai = 0.26 ± 0.11. Using the SCM Equa- 
tions (Equations 1 and 2) and the inferred values of V50, 
V50 and /so for SN 2003gd discussed, we determine the dis- 
tance in each band to be D(V) = 9.52 ± 4.14 Mpc and 
D(I) = 9.67 ± 3.82 Mpc. A weighted average yields a value 
of D = 9.6±2.8 Mpc. This error is statistical and comes from 
combining the uncertainties of each parameter in the SCM 
equations. Hamuy (2004b) finds a dispersion in the Hubble 
diagram of 0.3 m which translates to a precision of 15% for 
extragalactic distances. The distance estimate of M74 has a 
greater uncertainty as we have adopted large errors to ade- 
quately cover the assumptions we have made. 

We could also estimate the distance of SN 2003gd by 
using the Cepheid distance to SN 1999em (Leonard et al. 
2003, 11.71 ± 0.99 Mpc) directly. If we assume both super- 
novae have the same absolute magnitudes we can use the 
distance modulus for NGC 1637 and apply a correction for 
the difference in apparent magnitudes, shown in Equation 3, 
to calculate the distance modulus of SN 2003gd. We use the 
results given in Table 4 and the visual extinction of §3.2. 



03gd 



03gd 
03gd 



03 gd 



M 

J\®3gd 



- M 

j^OSgd 



03gd 
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m 99em + j^em + ^ 
^99em + Am + j^em _ 



03gd 



M 
_ j^)3gd 

^03gd 



(3) 



Using this equation and the VRI data we calculate fi(V) = 
30.05 ± 0.57, niR) = 30.15 ± 0.56 and n(I) = 30.16 ± 0.54. 
Together these give a weighted average of /1 = 30.12 ± 0.32 
or D = 10.6 ± 1.6 Mpc which is in agreement with our SCM 
distance estimate. 
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4.2 Brightest supergiant distance estimate using 
HST photometry 

This method uses the correlation between the average lu- 
minosity of the brightest supergiants and the host galaxy 
luminosity to estimate the distance. This average luminos- 
ity should be independent of the host galaxy luminosity for 
there to be no distance degeneracy. This method has pre- 
viously been used by Sohn & Davidge (1996) and Sharina, 
Karachentsev & Tikhonov (1996) to estimate the distance to 
M74, yielding similar results. Sohn & Davidge (1996) used 
photometry of supergiants in the disk of M74 and found 
H — 29.3, whereas Sharina et al. (1996) calculated the mean 
modulus of the group of M74 and its dwarf companions as 
H — 29.46. Both of these studies used ground based imaging 
to determine the magnitudes and colours of the brightest 
stars. The resolution of these images hinders the viability of 
the method (as discussed in §4.3), hence we have repeated 
this method with HST photometry of M74. 

We have used the method of Sohn & Davidge (1996) to 
determine the magnitude of the brightest supergiants. Their 
method firstly divides the supergiants into blue and red us- 
ing their (V — R) colours. The authors used (V — R) = 0.5 
which, taking into account the foreground reddening of M74, 
corresponds approximately to an F8 supergiant (Drilling 
& Landolt 2000, Table 15.7). The luminosity functions of 
blue and red supergiants then tell us the luminosity of the 
brightest supergiants, whilst removing contamination from 
foreground stars. This is accomplished by choosing the bin 
with statistically more supergiants than foreground stars 
and adopting this magnitude for the brightest supergiants. 
The number of foreground stars were calculated from field 
11 of Bahcall & Soneira (1981). In this work we have used an 
excess, of supergiants-to-foreground stars, of 2a to indicate 
a significant detection. 

The HST (WFPC2, see §2.3.2) photometry was firstly 
dereddened using the foreground extinction towards M74 
from Schlegel et al. (1998) using the NED interface. 2 The 
supergiants were divided into blue and red using their (B — I) 
colours and the intrinsic colour of Drilling & Landolt (2000, 
Table 15.7) for F8 type supergiants i.e. (B — I) = 1.28. 
We have assumed that supergiants with (B — I) < 1.28 
are blue supergiants (BSGs) and supergiants with B — I ^ 
1.28 are red supergiants (RSGs). The supergiants were then 
divided into luminosity bins and counted, where ubsg and 
nnsG are the number of blue and red supergiants within the 
bin, respectively. The resulting luminosity functions for blue 
and red supergiants are listed in Tables 10 and 11, and are 
plotted in Figures 16 and 17, in B and I respectively. 

The number of foreground stars were calculated from 
the predictions of Bahcall & Soneira (1981) and scaled to 
the WFPC2 field-of-view of 5.0 sq. arcmin. A significant sta- 
tistical excess of BSGs was found in the range B — 20 — 21 
and RSGs in the range I = 20 — 21. The B magnitude of the 
brightest BSGs was taken to be in the middle of this bin, 
B — 20. 5 ±0.5. Using the calibrations of Rozanski & Rowan- 
Robinson (1994, Figures 10(c) and 10(f)) for BSGs and the 
apparent B magnitude of M74 (Shostak & van der Kruit 
1984, B = 9.77), the distance moduli were calculated. Since 
their calibrations for RSGs are only for the K and V bands, 

2 http: / / nedwww.ipac.caltech.edu/forms / calculator.html 



Table 10. Luminosity functions of blue and red supergiants in 
B-band, where ngsG an d n RSG are the number of blue and red 
supergiants, respectively, and nf g is the number of foreground 
stars predicted from Bahcall & Soneira (1981), field 11. 



B 


n BSG 


nRSG 


n fg 


17-18 


1(1) 


0(0) 


0.3 


18-19 


0(0) 


1(1) 


0.5 


19-20 


0(0) 


0(0) 


0.6 


20-21 


6(2) 


3(2) 


0.8 


21-22 


71(8) 


10(3) 


0.9 


22-23 


336(18) 


37(6) 


1.1 


23-24 


1007(32) 


181(13) 


1.3 


24-25 


484(22) 


217(15) 


1.5 


25-26 


21(5) 


72(8) 


1.7 



Table 11. Luminosity functions of blue and red supergiants in 
/-band, where tlrsq and ursq are the number of blue and red 
supergiants, respectively, and nj g is the number of foreground 
stars predicted from Bahcall & Soneira (1981), field 11. 



I 


riBSG 


riRSG 


n fg 


17-18 


1(1) 


1(1) 


1.1 


18-19 


0(0) 


3(2) 


1.4 


19-20 


4(2) 


5(2) 


1.8 


20-21 


23(5) 


21(5) 


2.2 


21-22 


138(12) 


109(10) 


2.6 


22-23 


650(26) 


287(17) 


2.8 


23-24 


1060(33) 


95(10) 


2.8 


24-25 


43(7) 


6(2) 


2.8 


25-26 


7(3) 


0(0) 


2.8 



we converted our / magnitudes to V magnitudes. Rozanski 
& Rowan- Robinson warn against this and point out that the 
brightest stars in one band will not necessarily be the bright- 
est in another. We therefore chose the three bluest stars in 
the I — 20 — 21 bin. We averaged their B — I colours and 
their / magnitudes and used the intrinsic colours of Drilling 
& Landolt (2000, Table 15.7) to convert them to a V band 
magnitude. We obtain I = 20.6 with B — I = 1.3 which 



18 



20 



22 



24 



26 



Figure 16. Luminosity function in the B-band for BSGs where 
(B - I) < 1.28 in M74. 
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Figure 17. Luminosity function in the /-band for RSGs where 
(B - I) ^ 1.28 in M74. 



Table 12. Results of the brightest supergiant distance estimate. 



10(a) 


10(c) 


10(d) 


10(f) 


M74 29.70(0.92) 


29.47(1.03) 


29.40(0.90) 


29.15(1.01) 



corresponds to a F8-G0 supergiant with a (V — I) colour 
of 0.7 and magnitude of V = 21.4. We then used Figures 
10(a) and 10(d) of Rozanski & Rowan-Robinson (1994) to 
estimate the distance modulus for each calibration. The re- 
sulting distance moduli are listed in Table 12, where the 
errors are the combined errors from Rozanski & Rowan- 
Robinson (1994, Table 5), the error in the brightest super- 
giant magnitude and an error for the host galaxy reddening. 
We have also included an extra error of 0.5, to account for 
the transformation of the I magnitudes. A weighted aver- 
age gives n = 29.44 ± 0.48 corresponding to a distance of 
D = 7.7 ± 1.7 Mpc. Had we excluded the RSG result in light 
of the warnings of Rozanski & Rowan- Robinson (1994), the 
distance modulus would be /i = 29.31±0.72 or D = 7.3±2.4 
Mpc. The distance we find is consistent with those of Sohn 
& Davidge (1996, D = 7.2 Mpc) and Sharina et al. (1996, 
D = 7.8 Mpc). 

It is interesting however to note that Sharina et al. 
(1996) use the calibrations of Karachentsev & Tikhonov 
(1994) which are different to those of Rozanski & Rowan- 
Robinson (1994). Karachentsev & Tikhonov find three dif- 
ferent calibrations, two for BSGs and one for RSGs. These 
are shown in Equations 4—6, with a standard deviation of 
a(M) = 0.30, 0.30 and 0.45 respectively. 

(M v ) = 0.19M 3 B al -4.52 (4) 
(M B ) = 0.35M B al - 2.50 (5) 
(M B ) = -0.51(B - B 9ai ) - 4.14 (6) 

The distance modulus can be calculated from Equations 
4 and 5 using an iterative approach starting with the ab- 
solute magnitude from Shostak & van der Kruit (1984, 
M 9 B al = -20.23). Using this method we found the distance 
moduli to be n(V) = 29.71 ± 0.77 and n(B) = 30.12 ± 0.59, 
where the errors are the combined errors of the relation- 



Table 13. Summary of distance estimates. 



Method 


Source 


Distance 
(Mpc) 


Mean 


SCM 


1 


9.6(2.8) 


9.6(2.8) 


Brightest supergiants 


1 


7.7(1.7) 


8.1(2.3) 




1 


9.6(2.2) 






2 


7.2(2.0) 






3 


7.8(1.6) 




Kinematic 


4 


10.1 


10.2(2.6) 




5 


11.1 






6,7 


9.5 




Mean 






9.3(1.8) 



1 this paper, 2 Sohn & Davidge (1996), 3 Sharina et al. (1996), 
4 Tully (1988), 5 Ferguson et al. (1998), 6 Theureau et al. (1998), 
7 http://leda.univ-lyonl.fr 

ship, the error in the magnitude of the brightest supergiants 
and the reddening of the galaxy, which were both discussed 
previously. Sharina et al. (1996) used only the calibration 
given in Equation 6 for BSGs, which does not suffer from 
distance degeneracy. The authors found only one star red 
enough to be a RSG, but concluded it was a foreground 
star due to its magnitude. The distance modulus calculated 
from this relationship is /i(-B) = 30.11 ±0.67. A mean using 
the first two relationships only, gives a distance modulus of 
H — 29.91±0.49 corresponding to a distance of D = 9.6±2.2 
Mpc. This is a difference of 1.9 Mpc using the same data, 
but with two different calibrations. 



4.3 Summary and discussion of distance estimates 

A summary of distances to M74 from three different meth- 
ods is listed in Table 13. As well our SCM and brightest 
supergiant distance estimates, derived in §§4.1 and 4.2 re- 
spectively, we have collated three kinematic and two other 
brightest supergiant distance estimates from the literature. 
The kinematic distances shown in the table, use the heliocen- 
tric velocity of 656 km s _1 and apply three different models 
of infall onto the Virgo cluster. The values have been cor- 
rected for a Hubble constant of Ho = 72 km s _1 Mpc -1 . 
The simple mean of these distances is also listed in the table, 
where the error is a combination of the error on the mean 
and the uncertainty due to the cosmic thermal velocity dis- 
persion of 187 km s _1 (Tonry et al. 2000). The accuracy of 
these kinematic distances is limited not only by the observed 
velocity dispersion around the Hubble Flow, but on an ac- 
curate model to account for the infall onto Virgo and the 
Great Attractor. Depending on the model the values range 
from 9.5 — 11.1 Mpc, which is a 17% difference, hence there 
is a large uncertainty associated with this method. 

The brightest supergiants distances listed in Table 13 
come from three different sources. The first two distances, 
D = 7.7±1.7 and 9.6 ± 2.2 Mpc, are derived in §4.2. The 
third, D = 7.2 ± 2.0 Mpc, is the result of Sohn & Davidge 
(1996) where the error is a combination of the error in the 
size of the authors' bin and the error associated with the 
calibration used (Rozanski & Rowan-Robinson 1994, Figure 
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10a). The fourth, D = 7.8 ± 1.6 Mpc, is from Sharina et al. 
(1996) where the error is the combined error associated with 
the relationship used (Karachentsev & Tikhonov 1994) and 
the magnitude of the brightest BSGs. The simple mean is 
also listed in Table 13 where the error is a combination of the 
statistical error and the systematic error of Api = 0.59 which 
corresponds to a distance of 2.2 Mpc (Rozanski & Rowan- 
Robinson 1994; Karachentsev & Tikhonov 1994). Although 
the error in the average distance is the lowest of the three 
methods, the estimates range from 7.2 — 9.6 Mpc which is 
a 26% difference. There are large systematic uncertainties 
associated with this method and the derived distance also 
appears to be calibration dependent. 

Rozanski & Rowan-Robinson (1994) re-examine most 
of the existing observational data and investigate claims 
that the luminosity of the brightest BSGs exhibit a depen- 
dence on parent galaxy luminosity, introducing a distance 
degeneracy. Not only do they confirm this, a similar de- 
pendence is also found for the brightest RSGs, which was 
previously thought to be independent of galaxy luminos- 
ity. Rozanski & Rowan-Robinson also find that the result- 
ing error in the distance modulus from this method is far 
greater than previously believed. Whatever calibration they 
used, the error in the distance modulus was found never 
to be less than 0.55 mag. Piotto, Capaccioli & Bresolin 
(1992) and Karachentsev & Tikhonov (1994) also find such 
dependencies but find a smaller error in the regression of 
0.30 mag. Rozanski & Rowan- Robinson (1994) conclude that 
this method should not be used in determining the distance 
to single galaxies. Besides these calibration problems, this 
method also suffers from crowding, misidentification of su- 
pergiants with Hn regions, foreground stars, unresolved clus- 
ters and OB associations. These problems associated with 
poor resolution cause the distance to be systematically un- 
derestimated. This is illustrated by Leonard et al. (2003) 
who derive a Cepheid distance to NGC 1637 that is 50% 
larger than the brightest RSG distance estimate of Sohn & 
Davidge (1998). 

The estimates from this paper use HST photometry 
which has much better resolution than ground based obser- 
vations and should resolve individual stars from Hn regions 
and compact clusters, reducing this systematic error. If this 
systematic error is present we should find larger distances, 
which is what is found using the calibrations of Karachentsev 
& Tikhonov (1994). Although our first estimate is compara- 
ble to both Sharina et al. (1996) and Sohn & Davidge (1996) 
our second estimate using a different relationship differs by 
2 Mpc. It is worthy of note that whereas we have used all four 
of the calibrations of Rozanski & Rowan- Robinson (1994), 
Sohn & Davidge (1996) have used just that of Figure 10(a) 
for RSGs. If we look at the individual estimates of fi in Ta- 
ble 12 we see that fj,(a) = 29.70 is consistent with our results 
using Karachentsev & Tikhonov, whereas our other results, 
/i(c— >/), are not. If we use the parameters employed by Sha- 
rina et al. (1996) to calculate the distance modulus of M74, 
using the alternative relationships of Rozanski & Rowan- 
Robinson (1994, Figure 10c), we find fi = 28.92 instead of 
their \i — 29.3, which is a difference of 1.2 Mpc. Table 14 
shows the coefficients of the relationship M = aM 9 g l + b, 
found by each of the two different sources plus a third from 
Piotto et al. (1992). The coefficients of the relationship for 
RSGs, M(V), are consistent, but when we consider those of 



Table 14. Summary of coefficients from Piotto et al. (1992), 
Karachentsev & Tikhonov (1994) and Rozanski & Rowan- 
Robinson (1994) for M = aM g £ l + b. 



M a b o-(M) Source 



M(V) 


0.21 


-4.01 


>0.3 


1 




0.19 


-4.52 


0.30 


2 




0.21 


-4.10 


0.58 


3 


M(B) 


0.36 


-2.29 


~0.56 


1 




0.35 


-2.50 


0.30 


2 




0.28 


-3.45 


0.90 


3 



1 Piotto et al. (1992), 2 Karachentsev & Tikhonov (1994), 
3 Rozanski & Rowan-Robinson (1994) 



the BSGs relation, M(B), there is a disagreement between 
that of Rozanski & Rowan-Robinson (1994) and the other 
two papers. Rozanski & Rowan-Robinson find a much shal- 
lower relation resulting in the discrepancy between distance 
moduli. Given the agreement between the coefficients we 
consider the larger distance modulus to be more realistic. 

The SCM distance estimate also has its uncertainties. 
This method requires data at a phase of 50 days which is not 
always available. In order to obtain the necessary photomet- 
ric data, Hamuy (2004a) interpolated the observed V and 
I magnitudes to 50 days post explosion. In most cases the 
author was also able to interpolate the velocities by fitting 
a power-law to the observed velocities (Hamuy 2001), to an 
accuracy of ±300 km s~ . In four cases though it was nec- 
essary to extrapolate, introducing an error of 2000 km s _1 . 
The redshifts used in Hamuy (2004a) were derived from the 
observed heliocentric redshifts and converted to the Cos- 
mic Microwave Background frame in order to remove their 
peculiar motion. This in itself is limited by the accuracy 
of the model and by the cosmic thermal velocity. Hamuy 
noted that the least satisfactory aspect of the SCM is the 
dereddening method based on the comparison of supernova 
colour curves, which yields different results depending on 
what colour is used. It also depends on the accuracy of the 
reddening determination of the model supernova. Whilst not 
as accurate as distances derived from SN la, the SCM does 
give distances to within 15%. This method also gives com- 
parable values of Ho to those of Freedman et al. (2001). 

In §3.3 we highlighted the difficulties in measuring the 
minimum of the blue-shifted trough of the Fe II A5169 line 
at later epochs. We found this was not as much of an issue at 
earlier epochs where the velocity is needed for the SCM. A 
large uncertainty however was introduced in our application 
of the SCM, due to our lack of data at 50 days. Although this 
is the case we have been conservative in our error analysis 
and the SCM distance of D = 9.6 ± 2.8 Mpc is consistent 
with other distance estimates and associated errors. 

Each of the three different methods to measure distance 
have similar uncertainties and all have shortcomings. We 
have therefore taken an unweighted mean of the results of 
each method to give an absolute distance of D = 9.3 ± 
1.8 Mpc, where the error is the error on the mean. 
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5 ESTIMATE OF EJECTED 5b NI MASS 

The comparison of the light curve of SN 2003gd with that 
of SN 1999em in Figure 4, highlighted a major difference 
between the two supernovae. The lower tail luminosity was 
attributed to a smaller amount of nickel synthesised in the 
explosion as was also suggested by Van Dyk et al. (2003). 
Here we use three different methods to ascertain the mass of 
56 Ni produced. It was first of all estimated using the method 
of Hamuy (2003), secondly using a direct comparison with 
the light curve of SN 1987A and lastly using the new though 
unconfirmed method of Elmhamdi et al. (2003a). 



5.1 Nickel mass from bolometric luminosity of 
exponential tail 

Hamuy (2003) derives Mm from the bolometric luminosity 
of the exponential tail, assuming that all of the gamma rays 
due to 56 Co — > 56 Fe are fully thermalised. We first convert 
our U-band photometry in the tail to bolometric luminosi- 
ties using the formula given in Equation 1 of Hamuy (2003) 
which is given here in Equation 7. The bolometric correction 
is BC = 0.26 ±0.06 and the additive constant converts from 
Vega magnitudes to cgs units (Hamuy 2001, 2003). We have 
used the distance derived in §4.3. 



M m 



0.0" x f-^rl M 



log 



Mi, 



erg s _1 



-(V-A v + BC) + 5 log D - 8.14 



7.866 x 10~ 44 Lexp 



2.5 

(t - t(,)/(l + z) - T Ni 



(7) 



M (8) 



The nickel mass was then found using Equation 2 of Hamuy 
(2003) given here in Equation 8, where to is the explosion 
epoch, TNi = 6.1 days is the half-life of 56 Ni and rco = 111.26 
days is the half-life of 56 Co. Using this method we estimated 
Mm for each of the points in the tail and by taking a sim- 
ple average of these we found M Ni = 0.016^;^ M . This 
is approximately a third of that produced in SN 1999em, 
hence the lower tail luminosity. There is a large uncertainty 
associated with this value as the nickel mass calculation is 
strongly dependent on the explosion time, extinction and the 
distance, all of which are not well known. We have therefore 
calculated the range of nickel masses, for each point, from 
the extreme values for each parameter in Equations 7 and 8. 
We have taken an average of the maximum values and min- 
imum values and have estimated a more appropriate error 
from these. 



5.2 Nickel mass from a direct comparison to 
SN 1987A light curve 

The 56 Ni mass was also estimated from the difference in the 
UVOIR light curves of SN 2003gd and SN 1987A, assuming 
the same 7-ray deposition. A x 2 -fitting algorithm was used 
to shift the light curve of SN 1987A onto that of SN 2003gd 
to find the best fit. When constructing the light curve of 
SN 1987A a distance of 50 kpc was adopted. The difference 
in log luminosity was found to be log(L 87A /L) = 0.664 ± 
0.193. Equation 9 was then used to scale the nickel mass of 
SN 1987A, which was taken to be 0.075 M (e.g. Turatto 
et al. 1998), to estimate that of SN 2003gd. 



\L» 7A ) 



(9) 



Using this method we find M Ni = 0.016 ± 0.008 M where 
the error is the error associated with the fit combined with 
an error of 0.003 M related to the uncertainty in the ex- 
plosion date. This method is however still dependent on the 
distance to the supernova although the error in the nickel 
mass amply accounts for this. 



5.3 Nickel mass from "steepness of decline" 
correlation 

Elmhamdi et al. (2003a) reported a correlation between the 
rate of decline in the V-band, from the plateau to the tail, 
and the nickel mass estimated from the SN 1987A method. 
The authors defined a "steepness" parameter, S, which is 
the maximum gradient during the transition in mag d _1 . A 
sample of ten SNe II-P were used in determining the best 
linear fit which is given in Equation 10. 



log 



M Ni 



-6.2295S - 0.8147 



(10) 



The "steepness" parameter for SN 2003gd was found to be 
0.25 which yields the rather low result of Mm = 0.004l°;o! H ! 
Mq. The error was estimated from the scatter around the 
linear fit in the log(M N i/M )/5 plane, log(M Ni /M ) = 
0.147, and Equation 10 was used to translate this into a 
mass. This error is probably much larger due to the error in 
the fit. 



5.4 Discussion of nickel mass estimate 

The three independent methods, to obtain nickel mass, give 
two concurrent results of 0.016 M and one rather low es- 
timate of 0.004 Mq. The first two methods are both depen- 
dent on the explosion date and the distance, although the er- 
rors sufficiently accommodate these. The method of Hamuy 
(2003) assumes that the 7-rays are fully thermalised, but as 
the slope of the exponential tail of SN 2003gd follows the 
decay of 56 Co this is not unreasonable. The direct compari- 
son to the light curve of SN 1987A assumes that SN 2003gd 
deposited a similar fraction of 7-rays as SN 1987A. As the 
late-time decline of SN 1987A was very close to the decay 
rate of 56 Co, this assumption is also not unfounded. The 
method of Elmhamdi et al. (2003a) has the advantage of be- 
ing independent of explosion epoch and distance although 
it is an unconfirmed method. The correlation itself is dis- 
tance dependent and relies on the SN 1987A method. The 
result is very low and would suggest that the adopted dis- 
tance of 9.3 Mpc is out by a factor of two, which is un- 
likely. Given that the first two results are concurrent, and 
the errors adequately reflect the uncertainties, we take the 
average of the first two results and estimate a nickel mass of 
Mm = 0.016±g;8Jg M . 



6 MODEL COMPARISON 

We compared the data of SN 2003gd with the output of a 
semi-analytic light curve code (see Zampieri et al. 2003, for 
a detailed description of the code and the adopted initial 



© 2004 RAS, MNRAS 000, 1-23 



A Study of the Type II-P Supernova 2003gd in M74 17 



conditions). The goal of this comparison is twofold. First of 
all, we want to put some constraints on the input parame- 
ters (in particular the opacity) of the semi-analytic code by 
checking to see whether the mass estimate from the code is 
consistent with the observed measurement of the progenitor 
mass of 8±2 M Q (Smartt et al. 2004). Secondly, we would 
also like to use the code to derive some limits on the SN dis- 
tance and explosion date, which are rather uncertain. The 
code performs an estimation of the ejected envelope param- 
eters from a simultaneous comparison of the observed and 
computed light curve, photospheric gas velocity and contin- 
uum temperature. The main input parameters of the enve- 
lope are listed in Table 15. The fraction of the initial energy 
that goes into kinetic energy and the colour correction fac- 
tor (that measures the deviation of the continuum radiation 
temperature T c from the blackbody effective temperature 
T e ff) are input constants and are fixed at the values 0.5 
and 1.2, respectively. An important physical quantity that 
must be specified and that may affect in a significant way 
theoretical mass estimates is the envelope gas opacity n. 

The parameters of the model fits, corresponding to dif- 
ferent choices of the explosion epoch and distance, are re- 
ported in Table 15, where the estimated error in the pa- 
rameters is ~20-30%. For models A and B we assume that 
the phase identification derived in §3.1 is correct, implying 
that the explosion epoch is JD 2452717 (86 days before dis- 
covery). The two models differ in the adopted SN distance, 
D = 7.3 Mpc for model A and D = 9.2 Mpc for model B. 
Model C assumes JD 2452737 (~65 days before discovery) 
for the explosion epoch and D — 7.3 Mpc for the distance. 

The results for model A are shown in Figure 18. The 
overall agreement with the observed light curve, photo- 
spheric gas velocity and continuum temperature is rather 
good. The data are reproduced by a slightly under-energetic 
explosion (E = 1.4 x 10 51 erg) that has ejected M = 11 M Q 
and Mm — 9 x 10 -3 Mg. The comparison of the inferred 
value of M and the observed measurement of the progenitor 
mass shows that the star was probably in the upper end of 
the estimated mass range and sets an interesting constraint 
on the average opacity of the ejected envelope gas. In fact, 
values of k ^ 0.3 cm 2 g _1 result in ejected envelope masses 
in excess of 11 M@, barely consistent with the progenitor 
mass estimated from direct pre-explosion observations, once 
the compact remnant mass is taken into account (Smartt 
et al. 2004). Thus, large opacities are needed, in the range 
k = 0.3- 0.4 cm 2 g" 1 (k = 0.34 cm 2 g" 1 in Figure 18). 
Opacities in this range are appropriate for an envelope com- 
prised mostly of hydrogen, with a composition close to solar. 
This suggests that little mixing has occurred in SN 2003gd, 
the helium and other elements synthesised during nucleosyn- 
thesis remaining buried beneath the envelope. 

Model B tests the agreement between the model and 
the data assuming a different SN distance. This model has 
a rather typical explosion energy (E = 1.8 x 10 51 erg), while 
the ejected Ni mass (Mni = 1.5 x 10~ 2 Mq) is about 5 times 
lower than the "average" value for a normal type H-P SN 
(Mm = 7.5 x 10~ 2 Mq). We note, however, that there is 
evidence to suggest a continuous distribution of Mm in type 
II-P SNe (see Figure 23). The ejected envelope mass in this 
model is M — 13 Mq . Considering the mass of the compact 
remnant, this result indicates that M may be too large to 
be consistent with a progenitor of ^ 12 Mq. Therefore, if 




50 100 




50 100 

t (days) 

Figure 18. Results for model A showing V-band luminosity, L, 
velocity of metal (Sc II A 6246) lines, V g> i, and continuum tem- 
perature, T c , as a function of time for SN 2003gd. The adopted 
distance modulus and the estimated total reddening arc fj, = 29.32 
and A v = 0.42. The assumed explosion date is JD 2452717. The 
solid lines represent the curves from the semi-analytic model of 
Zampieri et al. (2003). 



the phase of the SN at discovery is ~90 days, the distance 
should be probably lower than 9 Mpc, in the lower end of 
the estimated range of distances. 

The uncertainty in the explosion epoch, however, is 
rather large and it cannot be ruled out that the discovery 
epoch is actually closer to the explosion date than estimated. 
Model C was run to test this possibility. We find that the 
data of SN 2003gd are reproduced by a fairly low energy ex- 
plosion (E = 0.7 x 10 51 erg) with the ejection of 6 M and 
a small amount of Ni (Mm — 7 x 10 -3 Mq). In this case 
smaller values of the opacity and/or larger distances would 
make the model still in agreement with the observations. 
However, an event of this type would have a plateau lasting 
^ 90 days, which is shorter than what is commonly observed 
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in type II-P SNe. Furthermore, the theoretical photospheric 
gas velocity appears to overestimate significantly (>20%) 
the observed line velocity at discovery. 



7 DISCUSSION 

7.1 SN 2003gd - a typical type II-P supernova? 

Unfortunately SN 2003gd was discovered ~86 days after ex- 
plosion (§3.1) when it appeared from behind the sun. Since 
valuable information about this supernova was lost we have 
made various assumptions about the supernova based on 
the striking spectroscopic similarities with SN 1999em and 
their comparable progenitor masses. Although SN 2003gd 
is spectroscopically very similar to SN 1999em it is fainter 
in its tail luminosity due to the lower nickel mass synthe- 
sised in the explosion. Hamuy (2003) presented evidence to 
suggest that the plateau luminosities and velocities are cor- 
related with the amount of nickel produced (see also Figure 
23). The lower tail luminosity may therefore suggest that 
SN 2003gd was intrinsically much fainter than SN 1999cm, 
although the spectral similarities and velocities at the end 
of the plateau phase would suggest otherwise. A number of 
papers have recently featured peculiar low-luminosity, low- 
nickel type II-P SNe (e.g. Zampieri et al. 2003; Pastorello 
et al. 2004a, b), but the question is, could SN 2003gd also 
belong to that group and if so, were our earlier assumptions 
justified? 

Elmhamdi et al. (2003a) and Pastorello et al. (2004b) 
noted a rapid excess in the (B — V) and (V — R) colours of 
the low- luminosity, or "faint", SN 1997D at the end of the 
photospheric phase (see Figure 19). This behaviour was also 
observed by Pastorello et al. (2004b) in SN 1999eu, another 
"faint" supernova, in the form of a sharp spike. Pastorello 
et al. suggested that the excess in colour may be charac- 
teristic of this low- luminosity sub-group. Figure 19 shows 
the intrinsic colour curves of SN 2003gd alongside those of 
the prototypical peculiar "faint" SN 1997D (Turatto et al. 
1998; Benetti et al. 2001; Pastorello et al. 2004a) and the 
model "normal" SN 1999em. The (B — V)o colour curve of 
SN 2003gd does rise more rapidly than that of SN 1997D, 
but it does not reach a significant excess. The SN 2003gd 
colours of Van Dyk et al. (2003) also do not display this 
steeper rise. The (V — R)o colour evolution of SN 2003gd is 
slightly more rapid than that of SN 1999em, at the end of the 
plateau phase, but it is nearer to the evolution of the "nor- 
mal" supernova. There may be a small spike in the colour 
of SN 2003gd at the end of the plateau phase at around 120 
days, but the photometry is sparse and it is therefore diffi- 
cult to confirm its existence. We would also note that there 
may also be a slight peak in the colour curve of SN 1999cm 
at ~8 days later. From these colour evolutions SN 2003gd 
appears to be more like a "normal" type II-P supernova. In 
order to further investigate our assumptions we compared 
the UVOIR light curves and spectra of SN 2003gd with ex- 
amples of "normal" and "faint" type II-P supernovae. 

The UVOIR light curve of SN 2003gd is shown in Fig- 
ure 20 along with those of three other type II-P supernovae, 
SNe 1999em, 1999gi and 1997D where SN 1999gi is a "nor- 
mal" plateau event. These UVOIR light curves were calcu- 
lated by summing up all the contributions from the U BVRI 



□ □ □ 




* 2003gd 

A 2003gd 2 

O 1999em 3 

□ 1997D 4 



50 100 150 200 

Time since explosion (days) 

Figure 19. Intrinsic (B — V) and (V — R) colour curves of 
SN 2003gd alongside those of the prototypical peculiar "faint" 
SN 1997D and the "normal" SN 1999em. The superscripts in the 
figure denote the source of the photometry: (1) this paper, (2) 
Van Dyk et al. (2003), (3) Hamuy et al. (2001) and (4) Benetti 
et al. (2001). 
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Figure 20. Comparison between the Ltjbvri light curves of 
SNe 2003gd, 1999em, 1999gi and 1997D which have nickel masses 
of 0.016, 0.048, 0.022 and 0.006 M Q , respectively. The nickel 
masses for SNe 1999em and 1999gi are from Hamuy (2003), and 
the nickel mass for 1997D is from Pastorello et al. (2004b), which 
are corrected for the revised distances used here. The superscripts 
in the figure denote the source of the photometry: (1) this paper, 
(2) Hamuy ct al. (2001), (3) Leonard et al. (2002b) and (4) Benetti 
et al. (2001). 
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Table 15. Theoretical parameters of SN 2003gd from the semi-analytic model of Zampieri et al. (2003) 



19 





Ro 


M 




Vo 


E 


ft 




T eff 




(10 13 cm) 


(M ) 


(M ) 


(10 8 cm s" 1 ) 


(10 51 erg) 


(cm 2 g- 1 ) 


(days) 


(K) 


model A 


2.3 


11 


9 x 10~ 3 


3.2 


1.4 


0.34 


45 


4600 


model B 


2.7 


13 


1.5 x 10~ 2 


3.4 


1.8 


0.34 


40 


5000 


model C 


2 


6 


7 x 10~ 3 


3.2 


0.7 


0.34 


45 


4600 



Ro initial radius of the ejected envelope 
M ejected envelope mass 
Mjm; ejected 56 Ni mass 



E initial thermal+kinetic energy of the ejecta 
k gas opacity 

£rec,o time when the envelope starts to recombinc 



Vo velocity of the envelope at the outer shell T e f f effective temperature during recombination 



Table 16. Distance estimates of NGC 3184 found in literature. 



Method 


Source 


Distance (Mpc) 


Mean 


Tully-Fisher 


1 


14.9 


11.1 




2 


7.2 




Tertiary indicators 


3 


7.9 


7.9 


EPM 


4 


11.1 


11.1 


Kinematic 


1 


10.6 


9.9 




5 


9.1 




Mean 






10.0 ±0.8 



1 http://leda.univ-lyonl.fr, 2 Pierce (1994), 3 dc Vaucoulcurs 
(1979), 4 Leonard ct al. (2002b), 5 Tully (1988) 



observations of each supernova. The few observations in U 
were used to estimate the contribution to the total flux from 
that band. In the case of SN 1999gi, where there was no U- 
band photometry, the U -band contribution was taken to be 
the same as that for SN 1999em. The distances used to cal- 
culate the UVOIR light curves were 11.4±0.9, 10.0±0.8 and 
15.3 ± 1.9 Mpc for SNe 1999em, 1999gi and 1997D respec- 
tively. The distance used for SN 2003gd is that discussed in 
§4.3. The distance of SN 1999em used here is the straight 
average of those discussed by Leonard et al. (2003, Table 
9). The distance of SN 1999gi, who's host galaxy was NGC 
3184, was taken to be the simple mean of the distance es- 
timates found in the literature, which are summarised in 
Table 16 where the error is the error on the mean. A liter- 
ature search for the host galaxy of SN 1997D, NGC 1536, 
returned two kinematic distances, 14.0 and 16.7 Mpc, from 
Tully (1988) and LEDA, which yield an average of 15. 3± 1.3 
Mpc. All kinematic distances were corrected for a Hubble 
Constant of Ho = 72 km s _1 Mpc -1 . The plateau luminos- 
ity of SN 2003gd is similar to that of SN 1999em at the end of 
the plateau, but dips to below both SNe 1999em and 1999gi 
during the tail phase. Although the decline from the plateau 
to the tail is greater, and steeper, in SN 2003gd than in the 
normal type II-P supernovae, the tail luminosity is much 
brighter than that of the faint, nickel-poor SN 1997D. 

The spectrum of SN 2003gd at 100 days is compared to 
those of SNe 1999em and 1997D at a similar phase, in Fig- 
ure 21. The spectra have been flux scaled using the method 
of §2.2, with the spectra of 1997D being scaled by an extra 
factor of five for clarity. The spectrum of SN 2003gd shows 
some similarities to both SNe 1997D and 1999em, but more 
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Figure 21. Comparison between the spectra of SNe 2003gd, 
1999em and 1997D which have nickel masses of 0.016, 0.048 
(Hamuy 2003) and 0.006 (Pastorello et al. 2004b) M , respec- 
tively, corrected for the distances used in this paper. SNe 1999cm 
and 1997D are representative of a "normal" and a "faint" type 
II-P supernova respectively. 



so to that of SN 1999em with its broad spectral features. 
The spectral flux of SN 1997D is much fainter than that of 
SN 2003gd, and also has much narrower spectral features in- 
dicative of the lower velocities. The appearance of the spec- 
tra, the higher velocities, the colour evolution and higher 
plateau luminosities of SN 2003gd suggest that it is more 
akin to the "normal" type II-P events. 

7.2 Implications for the progenitor of SN 2003gd 

Smartt et al. (2004) present the discovery of a red supergiant 
which exploded as SN 2003gd. The galaxy was imaged six to 
nine months before the supernova explosion by HST and the 
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Gemini Telescope. HST images after the explosion confirmed 
the positional coincidence of the supernova with a single 
resolved star. The observational properties of the star were 
compared with stellar evolutionary tracks. The progenitor 
star's colours and luminosity were found to be consistent 
with those of a red supergiant with an initial main sequence 
mass of 8I2 Mq. However, the reddening and distance used 
in Smartt et al. (2004) has since been updated to Av = 
0.43±0.19 (§3.2) and D = 9.3±1.8 (§4.3), and we recalculate 
the bolometric luminosity of the progenitor star using these 
values. We find Mboi = —6.0 ± 0.7 with no difference in 
the corresponding luminosity of log(L/L Q ) = 4.3 ± 0.3. We 
therefore find no changes to the authors' conclusions. 

We use Equation 2 of Nadyozhin (2003) , shown here in 
Equation 11, 

iog (£) = o - 234Mv+2 - 9iiog (d^i) 

+ ^(loiSrO- 1 - 829 < n > 

to compare the observed properties of SN 2003gd with the 
observed mass of 8la M . Since we can only make assump- 
tions about the length of the plateau, we have used this equa- 
tion to constrain the plateau length to ascertain whether it is 
consistent with other type II-P SNe. Based on the observed 
mass we have estimated the ejected mass to be 6 Mq , leaving 
2 Mq for the mass of the compact remnant plus other mass 
losses from the system. We find that the observed mass is 
consistent with a plateau length of 67^25 days, where the er- 
rors are the combined errors in each of the parameters in the 
equation. For the error in the ejected mass we have added, 
in quadrature, an extra error of 1 Mq to the error in the 
observed mass. Using the full range of the observed mass, 
and other parameters, we find the supernova parameters are 
consistent with a plateau range of 40-119 days, although re- 
alistically this is unlikely to be below around 90 days, from 
what is commonly observed for SNe II-P. Within the errors 
this mass appears not to be inconsistent with the observa- 
tions and may suggest a mass in the upper end of the mass 
range. 

The model of Zampieri et al. (2003), when applied to 
SN 2003gd (see §6), also seems to favour a main sequence 
mass in the upper end of the observed range, although this 
is achieved with a distance in the lower end of our distance 
range (Model A). Model B estimates an ejected mass of 
13 Mq using 9.2 Mpc, so is directly comparable with the 
observed mass. This is just outside the observed mass range, 
but is not necessarily inconsistent given the theoretical un- 
certainties associated with the model. The model does ap- 
pear to rule out an explosion date closer to discovery (Model 
C). We would like to stress though that more observations 
of progenitors are needed to fully test and constrain this 
model. 

The observations of progenitors have important impli- 
cations for the low-luminosity "faint" SNe II-P. There are 
two very different plausible progenitor models for these, one 
being the low-energy explosion of massive stars (e.g. Tu- 
ratto et al. 1998; Benetti et al. 2001; Zampieri et al. 2003). 
In this model the collapsing core forms a black hole and a 
significant amount of fall-back of material occurs. An alter- 
native scenario is the low-energy explosion of low-mass stars, 
presented by Chugai & Utrobin (2000) who successfully re- 



produced the observations of SN 1997D with an explosion 
of 10 50 ergs and an ejected mass of 6 Mq. The findings of 
Zampieri et al. (2004) support the high-mass progenitor sce- 
nario. The authors find a bimodal distribution of supernovae 
in the nickel mass, initial main sequence mass (Mms) plane. 
A reproduction of the authors' Figure 1 (right) is shown here 
in Figure 22 with the observational points of SNe 2003gd and 
1999em added with filled circles. Zampieri remarks that the 
theoretical estimates of Mms from Zampieri et al. (2004), 
shown with open circles in the figure, will most probably 
need to be revised using the value of the opacity from §6. In 
this picture the upper branch forms the "normal" type II-P 
branch and the lower the "faint" branch, with SN 2003gd 
occupying the lower end of the "normal" branch. The alter- 
native low-mass scenario positions the "faint" SNe at the 
lower end of the "normal" branch in one continuous distri- 
bution. Until there are direct detections of the progenitor 
stars of these supernovae, we cannot determine if either of 
these models is a true model. 

SNe II-P produce a large range of nickel masses from 
around 10~ 3 to 10 _1 M . Hamuy (2003) showed that a cor- 
relation exists between the amount of nickel produced in the 
explosion and the absolute magnitude in the U-band dur- 
ing the plateau phase. A reproduction of Hamuy's Figure 3 
is shown here in Figure 23 with SN 2003gd shown with a 
filled triangle. Four other supernovae, SNe 1994N, 1999eu, 
2001dc and 2003Z represented with open stars, are from Pa- 
storello (2003). These four supernovae plus SN 1999br are 
from the group of "faint" type II-P supernovae. In this fig- 
ure, SN 2003gd is in the centre of a continuous distribution 
of supernovae between the low nickel, faint SN 1999br and 
the high nickel, bright SN 1992am. 

Smartt et al. (2003a) made a comparison between all the 
known, direct information available on supernova progeni- 
tors with the theoretical predictions of Heger et al. (2003) , 
whose work on the pre-supernova evolution of massive stars 
produced supernova population diagrams as a function of 
mass and metallicity. These models do not consider binary 
stars or rotation. We plot this supernova population diagram 
in Figure 24 and add SN 2003gd, which has a solar metallic- 
ity (Smartt et al. 2004) . We repeat the caution which Smartt 
et al. issues, emphasising that at present this is a qualita- 
tive comparison. Whilst saying this, the progenitor detected 
for SN 2003gd is consistent with the theoretical II-P region 
undergoing a O/Ne/Mg core collapse albeit close to the the- 
oretical boundaries. 



8 CONCLUSIONS 

We have presented photometric and spectroscopic data on 
the type II-P SN 2003gd. We compared the BVRI light 
curves with those of the well observed, spectroscopically sim- 
ilar SN 1999em using a x 2 -fitting algorithm. This analysis al- 
lowed us to estimate an explosion epoch of JD 245271 7±21, 
corresponding to a date of 18th March 2003, and the VI 
magnitudes during the plateau phase. It was noted dur- 
ing this analysis that the tail luminosity of SN 2003gd was 
much lower than that of SN 1999em. Using this explosion 
epoch we were able to compare the velocity evolution of 
SN 2003gd with those of other similar type II-P SNe. The 
velocity evolution was found to be comparable at all epochs 
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Initial mass (solar masses) 

Figure 24. Theoretical SNe populations from the models of Heger et al. (2003). All objects below the dotted line should produce low 
luminosity SNe due to black hole formation and subsequent fall-back of material onto the iron core. All objects to the right of the solid 
black line should have lost their hydrogen envelope by the point of collapse resulting in hydrogen-free SNe. The white regions indicate 
areas where SN should not be observed. The observed positions of the SNe are approximately estimated and placed accordingly, where 
the arrows indicate those with limits. The zero point on the y-axis is assumed to be log(-Z/Z©) = — 1. 
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Figure 22. Reproduction of Figure 1 (right) of Zampicri ct al. 
(2004) showing the bimodal distribution of type II-P supernovae, 
in the Mn;/Mms plane, from their model. The open circles show 
theoretical data from the semi-empirical model of Zampicri et al. 
(2003), where Mms is estimated from the ejected envelope mass. 
The filled circles show observational data for SNe 2003gd and 
1999em, where Mn; is estimated from the light curve and M^g m 
is an upper mass limit only (Smartt et al. 2002a). SN 2003gd 
inhabits the lower end of the "normal" type II-P branch. 



in our dataset. This enabled us to extrapolate the velocity 
evolution of SN 2003gd backwards to find the velocity during 
the plateau phase. 

We calculated three distances to M74 using two dif- 
ferent methods, the Standardised Candle Method (SCM; 
Hamuy 2003, 2004a,b) and the Brightest Supergiant Method 
(BSM; Rozanski & Rowan-Robinson 1994; Karachentsev & 
Tikhonov 1994). Firstly, using the inferred velocity and VI 
magnitudes, we estimated a distance of D = 9.6 ±1.8 Mpc 



Figure 23. Reproduction of Figure 3 of Hamuy (2003) showing 
the correlation between the mass of 56 Ni synthesised in the super- 
nova with the absolute magnitude in the V-band at 50 days after 
explosion. The figure has additional points of SN 2003gd, shown 
with a filled triangle, and four "faint" SNe 1994N, 1999eu, 2001dc 
and 2003Z, shown with open stars, from Pastorello (2003). 



with the SCM. We then used HST (WFPC2) photometry 
to estimate the distances of 7.7 ± 1.7 Mpc and 9.6 ± 2.2 
Mpc using two different BSM calibrations. Using these three 
distances and other distances within the literature we esti- 
mated an overall distance of 9.3 ± 1.8 Mpc to M74. This dis- 
tance allowed the mass of 56 Ni to be estimated at 0.016 M©. 

A comparison was made between the observed data of 
SN 2003gd and the semi-analytical models of Zampieri et al. 
(2003). These models firstly tested the input parameters of 
the code and secondly ascertained whether the distance and 
explosion epoch were consistent with the observed progeni- 
tor mass. It was found that in order for the model to pro- 
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duce an ejected mass consistent with the observed progenitor 
mass the opacities needed were large, k = 0.3 — 04 cm -2 g -1 , 
suggesting little mixing occurred in SN 2003gd. The results 
of the models suggested that the progenitor was probably 
in the upper end of the progenitor mass range of 8i| 
and the explosion was not likely to have occurred much later 
than estimated in §3.1. 

The low tail luminosity of SN 2003gd sparked a de- 
bate about the nature of this supernova, whether it was a 
"normal" type II-P or whether it belonged to the group of 
"faint", low nickel SNe. To ascertain this we compared the 
colour evolution, pseudo-bolometric light curve and a spec- 
trum (100 days) of SN 2003gd with those of SNe belonging 
to the "normal" and "faint" groups. We also compared the 
observed mass of the supernova with theoretical models. Us- 
ing all the information available for SN 2003gd we conclude 
the following: 

(i) The colour evolution of SN 2003gd is similar to that 
of the "normal" type II-P SNe, although the reddening at 
the end of the plateau phase is slightly steeper it does not 
show a colour excess like the "faint" supernovae. 

(ii) The plateau luminosity of SN 2003gd is comparable 
to those of normal type II-P SNe although the decline from 
the plateau to the tail is greater. 

(iii) The tail luminosity of SN 2003gd, despite the larger 
decrease, is greater than the tail luminosity of the "faint" 
type II-P SNe. 

(iv) The spectra are much more similar to the "normal" 
type II-P SNe, with their broad spectral features, than the 
"faint" SNe, which are much fainter and have much narrower 
spectral features. 

(v) The velocities of SN 2003gd are comparable to those 
of the normal type II-P group. 

(vi) SN 2003gd lies in the centre of the correlated, contin- 
uous distribution of SNe II-P in the Mm/ My plane which 
has the low nickel, faint SN 1999br and high nickel, bright 
SN 1992am at its extremities. 

(vii) The observed progenitor mass of SN 2003gd lies at 
the lower end of the "normal" type II-P SN branch in the 
Mm/M-MS plane (Zampieri et al. 2003). 

(viii) In the SN population diagram of Heger et al. (2003) 
the observed progenitor mass of SN 2003gd lies in the type 
II-P region, which undergoes O/Ne/Mg core collapse, albeit 
close to the boundaries. 

We therefore conclude that the observations and observed 
properties of SN 2003gd are consistent with the observed 
progenitor star, a red supergiant of mass 8t.t M (Smartt 
et al. 2004), which resulted in a normal type II-P supernova 
at the end of its life. 
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